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Abstract 

For a heavy-light meson with a static heavy quark, we can explore the light 
quark distribution. The charge and matter radial distributions of these heavy- 
light mesons are measured on a 16 3 x 24 lattice at f3 = 5.7 and a hopping 
parameter corresponding to a light quark mass about that of the strange 
quark. Both distributions can be well fitted up to 4 lattice spacings (r ~ 0.7 
fm) with the exponential form wf(r), where Wi(r) = Aexp(— r/rj). For the 
charge(c) and matter(m) distributions r c ~ 0.32(2)fm and r m ~ 0.24(2)fm. 
We also discuss the normalisation of the total charge and matter integrated 
over all space, finding 1.30(5) and 0.4(1) respectively. 
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I. INTRODUCTION 



Lattice QCD has had considerable success in the understanding of the energies of few- 
quark systems. However, the spatial distributions of the quarks in these systems have 
received much less attention. The reasons for this are two-fold. Firstly, unlike energies, these 
distributions are not directly observable, but only arise in integrated forms such as sum-rules, 
form-factors, transition rates etc. Secondly, as will be seen later, their measurement on a 
lattice is more difficult and less accurate than the corresponding energies. In spite of this, it 
is of interest to extract lattice estimates of various spatial distributions and this is the aim 
of the present study. 

There have been earlier lattice measurements of radial distributions. However, they differ 
from the present work in several ways. For example, in Ref. M the authors are interested 
in the coupling between 5-mesons and the it, which involves a pseudo- vector coupling of 
the 7r to a single quark. In contrast, here it is the charge and matter distributions that are 
studied and these involve vector and scalar couplings. Scalar coupling to single quarks was 
also studied in Refs. [§], |J. However, there only the scalar sum-rule was evaluated, since 
that was primary a study of the dependence of the meson mass on the quark masses. 

A knowledge of spatial distributions can be utilized in a variety of ways. For example, 
the charge distribution of the light quark (q) in a heavy-light meson (Qq) can be used to 
check possible potential models, where the distributions are calculated from wavefunctions 
generated by some differential equation containing an interquark potential V. In such mod- 
els |4[] there are several uncertainties - the form of V, the form of the differential equation, 
how to include relativistic effects - some of which are tuned to ensure the model reproduces 
the correct (Qg)-energies. The latter can be in several partial waves Sin, P1/2, P3/2, ••• and 
can be either the observed energies of, for example, the 5-meson or the results of a lattice 
calculation. However, if also spatial distributions are known a priori from, say, a lattice cal- 
culation, then the uncertainties in such potential models will be reduced. Another way that 
the (Qq) charge distribution could be of use is in the understanding of multi-quark systems. 
In few-nucleon systems (e.g. 3 He, 4 He) it is found that the nucleon-nucleon correlations are, 
at short distances, very similar to that in the tu>o-nucleon system - with differences arising 
at large distances due to the different binding energies ||. This important observation can 
then be exploited in models of multi-nucleon systems, by assuming that the internucleon 
correlations are dominated by their two-nucleon counterparts. In the corresponding multi- 
quark case it would be of interest if a similar simplification were to arise. Therefore, as a 
first step in that direction, a knowledge of the two-quark radial correlation in the basic (Qq) 
case is needed - later to be compared with those in, say, the (Q 2 q) or (Q 2 q 2 ) systems. This 
probing of the (Qq) or (Q 2 q 2 ) structures could be carried to a more fundamental level by 
measuring the form of the underlying color fields that lead to the interquark potential V 
and the radial correlations. Such an extension would be analogous to the study of these 
fields in the (QQ) and (Q 2 Q 2 ) systems - see Refs. 

It should be added that the heavy-light system (Qq) is the quark model equivalent to 
the hydrogen atom (Pe~). Therefore, from a general point of view, it is of interest in any 
discussions comparing the properties of two-body sytems constructed from two particles one 
of which is very much heavier than the other. Also the interactions in the two cases have 
common features - the coulomb potential oc e 2 /r of the hydrogen atom versus the one-gluon 
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exchange oc a/r in the heavy-light meson. 

In section |T| the Maximal Variance Reduction method, crucial for extracting reliable 
results, is briefly discussed. In section [TTI| the formalism is introduced for evaluating the two- 
and three-point correlation functions C(2) and C(3). In section [IV] variational methods for 
analysing the lattice data are explained. In section [V] the results are given and in section |VT| 
some conclusions are made. 



II. MAXIMUM VARIANCE REDUCTION METHOD. 

It has been demonstrated that light-quark propagators can be constructed in an efficient 
way using the so-called Maximum Variance Reduction(MVR) method. Since this has been 
explained in detail elsewhere, for example in Ref. |7|, the emphasis here will be mainly on 
the differences that arise when estimating on a lattice the correlation functions C(2), (7(3) 
needed for measuring spatial charge and matter densities. In the MVR method the inverse 
of a positive definite matrix A is expressed in the form of a Monte Carlo integration 

Ajl = \\^ Mi e*p(-^* M\ (1) 

where the scalar fields are pseudo-fermions located on lattice sites For a given gauge 
configuration on this lattice, N independent samples of the fields can be constructed by 
Monte Carlo techniques resulting in a stochastic estimate of Ajf as an average of these N 
samples i.e. Ajl = (0*0-,). The N samples of the fields can be calculated separately and 
stored for use in any problem involving light quarks with the same gauge configurations. 

In QCD the matrix of interest is the Wilson-Dirac matrix Q — 1 — kM, which is not 
positive definite for those values of the hopping parameter k that are of interest. Therefore, 
we must deal with A = Q^Q, which is positive definite. Since M contains only nearest 
neighbour interactions, A - with at most next-to-nearest neighbour interactions - is still 
sufficiently local for effective updating schemes to be implemented. In this case the light- 
quark propagator from site i to site j is expressed as 

G q = Gfi = Qjl = ((QiM*<l>j) = (V#i>- (2) 

This is the key element in the following formalism. The Wilson-Dirac matrix also leads to 
an alternative form for the above light-quark propagator from site i to near site j 

G' q = = 75 < (QjMfi > 75 = 7s < Wi > 7s- (3) 

Later, it will be essential to use at some lattice sites operators that are purely local. This 
then restricts us to using at such sites only the fields that are located on single lattice 
sites. In contrast the ipi fields - defined as Qik4>k ~ are not purely local, since they contain 
fields on next-to-nearest neighbour sites. 

In the above, the term "Maximal Variance Reduction" comes from the technique applied 
to reduce the statistical noise in Eq. [|. The lattice is divided into two boxes (0 < t < T/2 
and T/2 < t < T) whose boundary is kept fixed. Variance of the pseudofermionic fields is 
then reduced by numerically solving the equation of motion inside each box. This allows 
the variance of propagators from one box to the other to be greatly reduced. However, in 
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the case of a three-point correlation two propagators are needed and this is best treated by 
choosing one of the points to be on the boundary of the boxes while the other two are inside 
their own boxes. Furthermore, the field at the boundary must be local to avoid the two 
propagators interfering with each other. This means that only the fields can be used on 
the boundary and there they couple to the charge or matter operator. For the points in the 
boxes, the temporal distance from the boundary should be approximately equal to give the 
propagators a similar degree of statistical variance. 



III. THE CORRELATION FUNCTIONS C(2) AND C(3,R) 

In this section an expression is given for evaluating the two-point correlation (7(2) - 
needed for extracting the basic Qq energies (m a ) and eigenfunctions(ff). This is then 
extended to the three-point correlation (7(3) in order to measure the radial correlations. 
Details are given in the Appendix. 



A. Two-point correlation functions C(2) 



Given the above light-quark propagators, then the two-point correlation functions (7(2) 
needed for extracting the energies of a heavy-light meson can be expressed in the four 
ways shown in Fig. |l]. These are all the same upto statistical errors, but their combination 
improves the overall measurement. In Figs. |I| a) and b), the heavy(static)-quark propagator 
from site (x, t) to site (x', t + T) is simply 



G Q (x, t ; x', t + T) = -(l + 74 )f/ Q (x,t,T)5 x , x ,, 



(4) 



where U Q (x,t,T) = l\J=o U±(x, t + i) is the gauge link product in the time direction. As 
shown in the Appendix this leads to the expression 



(7(2, T) = 2 (Re 
where d and e are Dirac spin indices. 



U Q V e *(x,t + T)0 e (x,t)+ £ 0*(x,t + T)^(x,£)] 

e=3,4 d=l,2 



(5) 



B. Three-point correlation functions C(3) 



The construction of the three-point correlation functions (7(3) follows that of (7(2) and 
are depicted in Fig. ||. Since a probe is now inserted at a distance r from the heavy-quark (Q), 
two light-quark propagators enter - one from Q to the probe and a second from the probe 
back to the Q. The purpose of the probe is to measure the charge or matter distribution 
at a definite point r. Therefore, only those light-quark propagators that involve the local 
basic field at r can be used, since the ip field contains contributions from fields at next- 
to-nearest neighbour sites and so is non-local. Later, when fuzzing is introduced similar 
restrictions will enter. In this work two probes are studied: i) 9(r) = 74 that measures the 
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charge distribution (actually the light quark charge normalised to 1) and ii) 0(r) = 1 that 
measures the matter distribution of the light quark. 

As shown in the Appendix, the same techniques introduced to evaluate C(2) can be 
extended to C(3) giving the overall three-point correlation function as 

C(3, -t 2 , t u r) = C Q (3) =F C*(3) = 2( X^Yg - XfY? >, (6) 

where the X q and can be expressed in terms of <j)(f) and ipUip respectively. In Eq.^| the 
relative sign enters for 0(r) = 74, since the q and q have opposite charges. 

IV. ANALYSIS 

The correlations of interest are essentially obtained from the ratio (C(3, T))/(C(2, T)) 
by projecting out the ground state expectation value. However, the latter is only achieved 
in the limit T — > 00. In practice, on a given lattice at the maximum possible values of T, 
the signal to noise ratio becomes large and effects from excited states are present. In order 
to reduce this contamination, a set of wavefunctions is constructed by fuzzing the original 
local wavefunction. These wavefunctions generate a better hadron operator where the Qq- 
meson is created and destroyed. Then, together with the original local form, they serve as 
a variational basis for analysing the data. 

A. The effect of fuzzing 

Fuzzing enters in two ways. 

1. Firstly, the basic links containing the gauge field have two fuzzings. In the stan- 
dard notation of, for example Ref. ||, Fuzzl has 2 iterations and Fuzz2 6 itera- 
tions. In both cases, the factor multiplying the basic link is f p = 2.5 

i.e. [A fuzzed link] = f p - [Straight link]+ [Sum of 4 spatial U-bends]. 

2. The pseudo-fermion field at a given lattice site r is considered to have three 
forms: 

(a) The basic form that is simply a function of the single lattice site r. 

(b) The field at r is an average of the fields on the neighbouring six lattice 
sites(i) i.e. 

i 

(c) The field at r is an average of the fields on the six lattice sites(j)that are 
two lattice spacings from r i.e. 

2 (r) = ^f/(F M ^2,r,r j )0(r j ). 

i 
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Therefore, only the basic field is local(L) - with <Pn=i and <Pn=2 being increasingly non- 
local. This means that in the calculation of the above three-point correlation functions, 
only the basic field should be used at r - the position of the probe insertion. There are 
now two reasons for this: 1) The field on the boundary at t = must be local. 2) The 
operator insertion must be local. This restriction does not occur elsewhere, so that the ip 
fields, which connect directly to the heavy quark, can be the fuzzed forms constructed from 
Qik4>N- This means that the two- and three- point correlation functions have the same size 
(3 x 3) of overall correlation matrix - LL, LF 1 , LF 2 , F\Fi, FiF 2 and F 2 F 2 . 

B. The variational method 

There are many ways of analysing the above correlations in order to extract the quantities 
of interest i.e. energies and wavef unctions. Here a variational method described in Ref. || 
is applied. 

Firstly, the two-point correlation data (7(2) are analysed to give the energies (m a ) and 
eigenvectors (v) for the Qg-system. These are then used in analysing the three-point corre- 
lation data C(3) to give the charge and matter densities. 

Consider the correlation function C(2, T) as an n * n matrix - upto 3*3 in this case with 
the elements LL, LFi, ....F 2 F 2 . Each element Cy(2, T) is then fitted with the form 

M 

C y (2, T) « (7^(2, T) = £ < exp(-m T)«?, (7) 

a=l 

where mi is the ground state energy of the heavy-light meson. The statistically independent 
matrix elements of C(2) are then fitted by varying the v"j and m a to minimise the difference 
between (7(2) and (7(2). 

We illustrate the procedure for the 2*2 case, where the (7(2) can be expressed as the 
product of three 2*2 matrices 

<7<2) = * (-P(-- r > exp( _° m2 T) )c, where c =(^«) (8) 

and c T is the transpose of c. In Ref. [^] the rows of c are called the v- vectors vf. Once the 
c matrix is known, any 2*2 correlation matrix C can be evaluated for the ground (excited) 
state wavefunction corresponding to the extracted eigenvalue mi{m 2 ) by reversing the above 
procedure to give 

C = (c T )~ 1 Cc- 1 i.e. C af , = u?C ij (T)u?, (9) 

where the v% are components of the u- vectors in Ref. |J. These u- vectors are the columns 
of the c _1 matrix and satisfy the condition 

< 4 = <W- (io) 

For the 2*2 case 

u 1 = [v 2 F1 , -v 2 L ]/det{c), u 2 = vl]/det(c). (11) 
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In the case where C is the above two-point correlation function C(2) and C in Eq. [7] is 
a perfect fit to C, then the operation in Eq. || would result in C being diagonal with the 
diagonal elements simply being exp(— m a a). Of course, in practice, the fit is never perfect 
and the off-diagonal elements of C are a measure of this goodness of fit. This will be 
demonstrated later. However, as pointed out in Ref. U, for other correlation functions 
there is no reason for C to be diagonal. 



V. RESULTS 

The results are presented in two distinct parts. Firstly, the two-point correlation function 
is analysed to give the ground and excited state energies and eigenvectors i.e. mi, m 2 in 
Eq. H and u 1 , u 2 in Eq. |9]. Secondly, these eigenvectors are used to extract the charge and 
matter radial distributions from the three-point correlation functions. 

The actual pure gauge configurations (20 in number) and the pseudofermion fields <j) (24 
per gauge configuration) were taken from the tabulation generated for the work of Ref. |12| . 



These are for a 16 3 x 24 lattice with /3 = 5.7 with the Sheikholeslami-Wohlert improved 
clover action with csw = 1-57 - corresponding to a lattice spacing of a « 0.17fm - and a 
hopping parameter k =0.14077. The latter corresponds roughly to the strange quark mass. 
This can be seen from ref |14| where the same parameters in the light-light system (qq) 
predict a vector meson to pseudoscalar meson mass ratio corresponding to strange quarks. 
More details can be found in Ref. M. 



A. Analysis of the two-point correlation function C(2,T) 

Essentially the energies (E), in lattice units, are extracted from the C(2,T) in Eq. |5| 
using 

£?[C7(2 I T)] = -ln[^|SL] asT^oo. (12) 

In Fig. |3|a) the basic C(2,T) are plotted for the three diagonal matrix elements where the 
fields are: i) purely local (L), ii) all Fuzzl, iii) all Fuzz2 - in the notation of subsection |1V A. 



They are all seen to be well determined with the errors only being significant for Fuzz2 with 
T > 10. In Fig. |3]b) the results E[Ca(2,T)] from Eq. |T2| are plotted for the three sets of 
diagonal matrix elements with i = L, Fl, F2. There it is seen that only E[Cfifi{2,T)] 
shows a clean plateau that extends from T = 5 to 9 with a value about 0.88(1). 

To combine these results by a variational calculation using Eq. [?| two numbers are fixed: 

i) M - the number of energies being included. Here, this is taken to be the same as the 
number of paths for each energy and results in the correlation matrices being square. 

ii) T min - the minimum value of T used in the fit. 

Here we consider four possibilities to check the dependence of the final results on this 
fitting procedure: 

Case 1): In Eq. [7|, the Cij(T) are defined in terms of 3 paths (i.e. i,j = 1,2,3) and 3 
exponentials (i.e. M — 3) with T min = 3. This includes the Local and both Fuzzl and Fuzz2 
paths. 
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Case 2): Same as Case 1) but with 2 paths and 2 exponentials. This includes only the 
Local and Fuzzl paths. 

Cases 3,4) are the same as Cases 1,2) but with T min = 4. 

Minimising the difference between the C(2,T) and £7(2, T) in Eq. ^ gives the parameters 
in Table |. These are surprisingly good fits, when it is realised that the errors on the C(2, T) 
are, in most cases, much less than 1%. However, only Case 3 gives xV n cfo/(2) < 1 and so 
this is the one that will be used in most of this study. 

In Fig. |^c) the results for E[C aa ) from Eq. [12] are plotted for the two sets of diagonal 
matrix elements Cn(2,T) and 622(2, T) for Cases 3 and 4. As expected, E[Cu(2,T)] « 
0.86(2) and .£[622(2, T)] ~ 1.24(5) - energies that are consistent with the values of m\ and 
rri2 in Table Q. As a check on the off-diagonal matrix element 6*12(2, T), we evaluate 

This is seen to be at the 1% level. The conclusion to be drawn from Fig. |3]c) is that 6 is, 
indeed, approximately diagonal with C a/ g(2, T) ~ exp(— m a T)5 a p. These results will serve 
as a comparison when analysing C(3, T, r) later. 



B. Analysis of the three-point correlation function for sum rules 

The charge and matter radial distributions F(r,Q) of the light quark in the Qq system 
are extracted using 

where = 74 or 1. However, before showing these radial distributions, it is of interest to 
first study the corresponding sum rules 

cisumiywQ _ Er(6(3, @, T, T)) _ (C SU7n (3, Q, T)) 

b Jc{2~TY) " (6(2, T)> ' (15) 

where J2r represents the sum over all spacial lattice sites at time t = - the time when 
the probe acts. The summation can be easily carried out exactly on the lattice. For the 
charge distribution this sum rule should, in the continuum limit, simply yield the charge of 
the light quark, which we have chosen to normalise to unity. For the matter distribution 
the situation is less clear - see Refs. |J, f|. 

In Fig. |a) we show F sum [C u (^,T)] for the three cases % = L, Fl, F2. There it is 
seen that, at large T, the sum rule for i = L becomes 1.05(5) and so is consistent with 
unity. However, the results for i = Fl, F2 appear to be somewhat greater than unity 
for T > 5 converging to 1.25(5) and 1.4(1). When these results are combined using the 
u- vectors in Table [I], then - as seen in Fig. f|(b) - the ground state sums F sum [Cn{l4:,T)] 
for both Cases 3 and 4 are dominated by the F sum [CFiFi(j4,T)]. Consequently, they are 
again distinctly greater than unity with both Cases 3 and 4 tending to about 1.30(5). This 
seems unavoidable since, as seen from Table p], the u vector components u Fl and u -p<2 are 
an order of magnitude larger than u\. In Fig. |]b) the sum-rule for the first excited state 
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F sum [C22{j4:,T)] is also shown. This appears to be approaching unity at about T ^7 - 8. 
However, the signal is swamped by the error bars at larger T. 

The correlations C^" m (3) are not so well diagonalised as the CV,(2) were forced to be 
earlier. A measure of this is plotted in Fig. |]b) as 

*-[C a (e,n 1 = fg||g . (is) 

This is seen to be at the 10% level - an order of magnitude larger than the corresponding 
R[Ci 2 {2, T)] in Eq. |T^. Therefore, it is seen that the 74 sum rule calculated in the above 
way has two undesirable features: 

1. The ground state sum rule F sum [611(74, T)] is not consistent with unity for large 
T being more like 1.30(5). Only in the continuum limit should the sum be unity. 
To some extent this must be expected, since in the present non-continuum situa- 
tion the lattice vector current is not conserved. In principle this can be corrected 
in various ways. Unfortunately, at the present value of (3 = 5.7, perturbative ex- 
pressions would exhibit poor convergenge and so be unreliable. Non-perturbative 
corrections are discussed in the recent review in Ref. [|Kj. There an improved 
vector operator is introduced and calculations performed in the quenched ap- 
proximation at /3 — 6.2 and 6.0. However, they find that, whereas the (3 =6.2 
results are satisfactory, those at (3 = 6.0 are not - suggesting that 0(a 2 ) discreti- 
sation errors are not small at the larger lattice spacing. The situation would be 
even worse at the value of (3 — 5.7 used here. Even so, it is illustrative to see the 
results at the higher /3's, since they may indicate what could possibly be expected 
at lower /3's. In Ref. |L0|, when the expectation value of the renormalised vector 



current is expressed as 

(J R ) = F v ( l4 (x)} } (17) 

where (74(2)) is the quantity evaluated in Eq. [15|, they find that F v « 0.8. 
Similar reductions are found in Ref. |l| for the axial vector operator. Of course, 
the above argument could be reversed to say that the present charge sum-rule 
measurement of 1.30(5) gives a non-perturbative estimate of F v as 0.77(3). Later, 



when individual radial contributions are extracted as in Eq. IT1, this continuum 



effect could be incorporated approximately by the renormalisation 

F[CU^T,r)]^ Fsum[cu ^ T)y (18) 

2. The off-diagonal terms, as illustrated by -R sum [612(74, T)} in Eq. [16], are not 
zero. This is potentially more disturbing than point 1), since we wish to identify 
^[611(74, T, r)] as the ground state charge distribution, and so any off-diagonal 
effects would be difficult to interpret. Also they cannot be easily renormalised 



away as in Eq. IS . 
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So far the basic wavef unctions uf have been determined, from Eq. [j], via the vf i.e. by 
considering only the two-point correlations C(2). In an attempt to overcome the two above 
problems with F sum [C( r )4 L) T)], the analysis of the data is now carried out not only by fitting 
C{2) but also some features of C(3). 

1. In analogy to Eq. [7] the sum rule data are fitted with the expression 

M M 

C™ m (3, 7 4,T) « (™(3, 7 4,T) = E E <exp[-m a t 1 ]x^exp[-m /9 (T - *0]«f, 

(19) 

where the u™ are from Table | - the earlier result of minimizing the CV, (2) - but 
the x a P are treated as free parameters. However, a restriction must be made on 
the values of T used in Eq. [19]. From Table [I] it is seen that T min must be at least 
4 to ensure x 2 (2)/n<fo/(2) is comparable to unity. Therefore, in C sum (3) each of 
the two propagators should be at least 4 euclidean time steps i.e. in Eq. [19] we 
should have ti > 4 and (T — ti) > 4. This means that only the C(3,T) data 
with T min {3) > 8 should be included in any fitting procedure. However, in the 
following a series of T min {3) values, ranging from 4 to 9, are used to see how the 
final results depend on T min (3). Therefore, in Case 3, i.e. 3*3 with T min {3) = 4 [8], 
this involves fitting 42 [18] pieces of data with 6 parameters and in Case 4, i.e. 2*2 
with T m j„(3) = 4 [8], this involves fitting 21 [9] pieces of data with 3 parameters. 
In Table [T| the results for the x 2 / n dof are shown separately for the fits to C(2) 
and C sum (3). There for Case 3 it is seen that - as with x 2 (2)/n<fo/(2) - the 
X 2 {3)/ndof{3) are also ~ 1 provided T min {3) > 6. This result can be compared 
directly with Fig. |]b), since it is simply an alternative analysis under similar 
constraints. It is seen that the values of x 11 = 1.35(5) and x 22 m 1 are indeed 
consistent with the asymptotic values of F Sum in Fig. |]b). On the other hand, 
x 12 cannot be compared directly with i? sum [6*12(74, T)\ in Eq. All that can be 
said there is that both analyses result in non-zero and negative values for x 12 . 

2. In an attempt to overcome this last point, in Cases 3' and 4', all three off diagonal 
terms x a/3 are fixed at zero. For a given T TO j n (3), this decreases x 11 but at the ex- 
pense of increasing x 2 (3)/^do/(3). The outcome is that for x 2 (3)/ n dof(3) < 1 with 
x 12 = and T min {3) = 8 we get x 11 is 1.29(3) - a number consistent with 1.30(5) 
from Fig. |]b). However, x 22 changes violently, even though the plot corresponding 
to Fig. go) is very similar - with, in particular, F sum [6*22(74, T)] again approach- 
ing unity and not zero as would be expected from the value for x 22 in Table [TT]. To 
check this last unexpected result, the analysis program was run using the model 
results C , ^ nm (3,7 4 ,T) instead of the lattice data C^ m (3, 74, T). The plot corre- 
sponding to Fig. f|b) now has for the excited state F sum [6*22(74, T)} m 0.005 for all 
values of T - consistent with Table [TI]. This difference with x 22 demonstrates the 
need to try to improve the x 2 (3)/nrf /(3) at smaller values of T min (3). It appears 
that the results for the excited states are very dependent on the v-vectors, since 
they involve delicate cancellations. Of course, the main concern is the consistency 
between x 11 and F sum [C*n(74, T)\ and this emerges intact. 
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3. In the above, the data for (7^(2, T) and (7j?" m (3, 74, T) were fitted separately. 
Therefore, now a combined fit is made using Eqs. |7] and [l^ i.e. for Case 3(3') with 
T m in = 4, T min {3) = 8, this involves fitting 54+18 pieces of data for (7(2) + (7(3) 
with 12+6(3) parameters. However, this has a completely negligible effect e.g. 
in Case 3, T min = 4, T min (3) = 8 the values of x 2 (2)/ra do/ (2), x 2 {3)/n dof (3) and 
X 2 (2 + 3)/n rf0/ (2 + 3) change from 0.647, 0.718, 0.673 to 0.652, 0.699, 0.669 and 
x 11 from 1.3274(249) to 1.3274(255). The conclusion is that for the sum rules 
there is no benefit in fine tuning the results by simultaneously fitting C(2) and 
C(3). 



C. Analysis of the three-point correlation function for radial distributions 

In Figs. |5]a) and |5]b) are shown, for r/a = 0, . . . , 5, the ratios 

^ Tir)] .(^ga. (20) 

based on Eq. |TJ] using the v- vectors of Case 3. These all exhibit, to a certain extent, a 
plateau as T becomes large. However, for r/a > 5 the error bars become very large beyond 
T = 8. No attempt will be made at this stage to extract the asymptotic values giving the 
radial distributions. 

The second procedure for analysing the radial distribution three-point correlation func- 
tions C(3,9,T,r) is similar to the one followed above but using the expression 

M M 

Cij(3, 74, T, r) w Cy(3, 74, T, r) = £ £ v?exp[-m Q ti]x^(r) exp[-m^(T - t0]^, (21) 

a=l /3=1 

analogous to Eq. [B|for the sum- rule but with all the r-dependence being put into the x a @(r). 
Two types of fit are made: 

1. The v- vectors, obtained by minimising the C(2), are used in Eq. |2l] and for each 
value r the x af3 (r) are varied to ensure a good fit to Cy(3, 74, T, r) by the model 
form Cij(3, 74, T,r). 

2. A simultaneous fit of Cij(2,T) and CV,(3, 74, T, r) using (7y(2, T) and 

(3, 74, T, r) of Eqs. [7| and || 

To achieve a % 2 (2 + 3)/ndo/(2 + 3) ~ 1 is now more difficult than the earlier C(2), C sum (3) 
fit, since the CV,(3, 74, T, r) have error bars that are smaller than the corresponding sum 
rule correlations. For example, the most important correlations i = j = Fl have the 
values (7(2, T = 4) = 0.0752(3) and (7(3, T = 8, r/a = 2) = 0.0000528(9) compared with 
C sum (3,T = 8) = 0.029(1) i.e. the errors on C(2) are w 0.5%, those on C sum {3) are w 5%, 
but those on (7(3, r) are only 2%. This is seen in Table fTJ, where in particular for r/a = 1 
the values of x 2 {3)/n ( iof{3) are all greater than unity. In this table it is also seen that 
the radial distribution of the ground state charge density, x 11 (r) , is well determined and is 
consistent with the plateaux in Figs. |5]a) and b). When the x ll (r) are plotted on a semi- log 
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scale as in Fig. |, the distribution for r /a < 4 is approximately a straight line suggesting 
that u> c (r) = A c exp(— r/r c ), where w 2 (r) = x u (r). The function w c (r) is then interpreted as 
a radial wavefunction. However, the other diagonal matrix element, x 22 (r), is much less well 
determined - as is seen in the lower half of Table [II. In fact, for some values of r, it appears 
to be slightly negative. But the actual values are very small and could well be consistent with 
zero. Only for Case 3 at r/a = is there a definite signal with x 22 (0) ~ 0.19. This suggests 
x 22 (r) is very sharply peaked compared with x u (r). However, it must be remembered that 
in Case 3 the off-diagonal terms x a @(r) are not forced to be zero and so the interpretation 
of x 22 (r) is not clear. In Case 3' where the off-diagonal terms x a/3 (r) are forced to be zero, 
no signal can be extracted at r = 0. There is certainly no sign of a node that would be 
expected of an excited s-wave. 

Fig. H also suggests that for r/a > 4 the function w(r) drops off faster than the above 
simple exponential. Such an effect is expected at sufficiently large r when the linear rising 
confining potential becomes important. Then as r — > oo the wavefunction should have an 
exp[— (r/rx) 3 / 2 ] asymptotic form. Unfortunately, for r/a > 5 the errors on the data become 
too large. So that the above observation rests completely on the r/a = 5 data, which itself 
is rather poor. Even so, there the errors are still sufficiently small to support this. It is 
planned to extend the present calculation to off-axis points. The r/a = 5 results can then 
be checked by performing simulations at x/a — 3, y/a = 4. The resulting data should be 
more accurate since each (x, ?/)-pair arises in 24 different ways -not just 6 as for the on-axis 
points with r ^ 0. 

The actual parameters defining w c (r) can be extracted in a variety of ways depending 
on T min (3) and the range of r values used in the fit. But as seen in Table [TV] they appear 
to be quite stable for Case 3 and 3' separately. For Case 3, A c pa 0.23(1) and r c /a pa 2.1(1) 
and for Case 3', A c pa 0.26(1) and r c /a ~ 1.9(1). Given these numbers, then estimates of 
the sum-rule are I c = / °° drw 2 (r) = TcA 2 r^. These are shown in the last column of Table [IV 
and are all 1.6(1) for Case 3 and 1.5(1) for Case 3'. The corresponding numbers for the 
sum rule obtained directly from the lattice and shown in Table || are 1.35(5) and 1.25(5) 
respectively. This difference between the two methods is not surprising, since the integrand 
in I c = Jq 00 drw 2 (r) is maximum at r/a pa 2 and half of the sum I c is lying inside this value 
of r. Therefore, lattice artifacts, present at small values of r, may play a role. These are 
expected to affect I c more than the direct results of Table (Tl - a point that can be checked 
by estimating x 11 away from axes e.g. not only at (±1,0,0), (0, ±1,0), (0,0, ±1) but also 
at (±1,±1,0) etc. 

In addition to fitting the x 11 ~ w 2 (r) with simply the two parameter function 
w c (r) = A c exp(— r/r c ), attempts were made with the three parameter function w c (r) = 
A c exp[—(r/r c ) p ). As seen in Table |Iy], this results in values of p that are somewhat greater 
than unity - as expected from Fig. ^ where the r/a = 5 point appears to drop below the 
earlier simple exponential. This also has the effect of decreasing the value of J c = J °° drw 2 (r) 
from the simple exponential value of irA 2 c r^,. However, since about 90% of I c comes from r 
values less than 5, this in practice can have only a minor effect. 

Assuming the simple exponential form, an estimate of the mean-square charge radius is 
(r 2 ) = 3r 2 in lattice units of a ~ 0.17fm. For Cases 3 and 3' this gives 0.38(3) and 0.32(1) 
fm 2 respectively. 
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D. Analysis of the three-point correlation function for the Matter radial distributions 



The previous two subsections have dealt with the charge radial distribution, where the 
operator 6 in Eq. ^ is 74 for probing the q charge. In this subsection = 1, which probes 
the q matter. 

In Table [V] the values of the x af3 in Cases 3 and 3' are given. Here it is seen that x 11 the 
matter sum rule for the ground state is more erratic than its charge counterpart in Table [II 

- a point clearly seen in the corresponding plot in Fig. |7|. A reasonable estimate for x 11 is 
0.4(1). This value is consistent with the corresponding estimates in Refs. and |7[] and for 
12 3 x 24 lattices. These were made by employing data from different hopping parameters 
(k) and using the identity 

x n = ^h) ? (22) 

where ami is the ground state energy and k the hopping parameter - see Ref. @. When the 
mi's correspond to the cases where the light meson is of about one and two strange quark 
masses, Refs. and [[7] give 0.34(8) and 0.31(6) respectively - consistent with the present 
value of 0.4(1). 

These larger values are also consistent with the following simple estimate - again using 
the above identity in Eq. |22|. If the Qg-meson mass (ami) is taken to be simply the sum of 
the quark masses arriQ + am q and k^ 1 = 8 + 2am q , then 

u = d(am Q + am) = 

d(8 + 2am q ) v ; 

Another reason for expecting x 11 < 1 also follows from a potential approach using the 
Dirac equation. If the q is treated as a particle in a potential generated by the Q, then its 
wavefunction will be of the form (f,g), where f(g) are the large(small) components of the 
Dirac wavefunction. The charge of the q will then be simply Iq = J dr[f 2 (r) + g 2 (r)], which 
by the normalisation will be unity. However, when the charge operator(7 4 ) is replaced by 
the matter operator (unity), then the correponding integral is now Im = J dr[f 2 (r) — g 2 (r)]. 
In the non-relativistic limit Iq = Im- But as relativistic effects enter i.e. g 2 increases from 
zero, then Im becomes less than Iq i.e. less than unity. An indication of this is seen in the 
lattice results in Table 10.3 in Ref. 0, where x 11 decreases from 0.21(8) to 0.11(5) as the q 
mass decreases from about two to one strange quark masses. This also shows - as expected 

- that the interquark potential is more than one-gluon-exchange (— a/r), since the latter 
results in g/ f = (1 — 7) /a, where 7 = yl — a 2 - a ratio that is independent of the light 



quark mass - see, for example, [13]. For a = 0.3 this gives g/f = 0.15 - a number that 
is much smaller than the ~ 0.75 suggested by the charge and matter sum-rules measured 
above on the lattice. 

In Table |VI] the matter densities are given in analogy with the charge densities of Ta- 
ble [Til]. Comparison of these two tables shows that the ground state matter distribution 
x (matter) = x^ decays faster than the corresponding ground state charge distribution 
x(charge) = x]}. At r = the two are comparable but, by r/a = 4, x(matter) is only 25% 
of x(charge). This is also seen in Figs. ^. Now the r/a = 5 data is even more uncertain 
than the earlier charge data. It is, therefore, not quoted. 
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Table |VII| shows the parameters A m and r m when w m (r) is parametrized as w m {r) = 
A m exp(— r / r m ) . There it is seen that A c , for the charge distribution, and A m are comparable, 
but that r m pa 1.55(5) compared with the charge range of r c /a ~ 2.0. This difference between 
r c and r m means the interpretation that w c (r) and w m (r) are both q wavefunctions is not 
so direct, since in the most naive models one would expect w c (r) = w m (r). 



VI. CONCLUSION 

In this paper a first study has been made of the radial structure of a heavy-light meson 
- the quark equivalent of the hydrogen atom. This can be considered as a partial extension 
of Ref. in which only the energies of heavy-light mesons were measured on a lattice. Here 
the emphasis is on the charge and matter radial distributions. So far these distributions 
have only been extracted for the ground state, with the extension to other partial waves - 
as in Ref. M - only now beginning to be studied. A further extension would be a study of 
the form of the underlying color fields that lead to these radial correlations. This would be 
analogous to the studies in Refs. || for the (QQ) and (Q 2 Q 2 ) systems. 

The main result is in Figure |6|, where it is seen that both the charge and matter radial 
distributions fall off more or less exponentially as A 2 exp[— 2r/r^. The amplitudes Aj are 
about the same with A c pa 0.26(1) and A m pa 0.29(1), whereas the charge range r c /a pa 1.9(1) 
is considerably longer than the matter range r m /a ~ 1-4(1). This difference is reflected in the 
spatial sum-rule, which is about 1.30(5) for the charge and 0.4(1) for the matter. The fact 
that the charge sum-rule is not unity, as would be expected from vector current conservation, 



can be attributed to the finite spacing of the lattice. As shown in Refs. |TJ] and [IT5| this can 
easily lead to 10-20% reductions. 

It should be added that there are other definitions of Qg-wavefunctions. Here the rel- 
ative wavefunction [w(r)] is assumed to be real with w 2 {r) giving the charge density - the 
quantity actually measured from (C(3,T)/C(2,T)). Another form can be extracted di- 
rectly from a two-point correlation function in which the operators at the sink and source 
are of spatial size r% and r 2 respectively. The ground state correlation can then be fit- 
ted with WBs( r i) w Bs( r 2) exp(— m\T) to give a WBs( r ), which can be identified with the 
Bethe-Salpeter wavefunction - see Ref. [[?]]. 

The above charge and matter radial densities are related by a simple fourier transform 
to the momentum space vector and scalar form factors [F VtS (q 2 )] of the 5-meson. Unfortu- 
nately, the present densities are not accurate enough over a sufficiently large range of r to 
perform this transform. However, a simple model of these form factors is that they are dom- 
inated by the pole due to the lightest meson of mass m i.e. they have a form oc (q 2 + m 2 )~ 1 . 
On the lattice with our parameters and quark mass, it is found that the lightest vector and 
scalar meson masses are am v = 0.815(5) and am s = 1.39(5) respectively - see Ref. [I4|j9"|. 
The fourier transform of these pole terms is then a Yukawa form oc exp(— mr)/r, which - in 
principle - can now be compared directly with the charge and matter radial densities mea- 
sured here. However, it is only the asymptotic form that should be used in this comparison 
since that will be controlled by the lightest particle mass. But from Fig. |^ it is seen that the 
present data only extends up to r/a = 5(4) for the charge (matter) density, corresponding 
to mr = 0.4(0.6) respectively which is not large. Furthermore, the charge density is already 
well described by a simple exponential up to r/a = 4 and so a comparably good fit over this 
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range by a Yukawa form is ruled out. As a compromise, since the charge density data is 
somewhat better than that of the matter density, the charge densities at r /a = 3 and 4 are 
used to extract estimates of m v and the matter densities at r/a = 2 and 3 for an estimate 
of m s . Case 3 leads to am v =0.7(l) and am s = 1.4(4) - the corresponding numbers for Case 
3' being 0.8(1) and 1.1(1). Even though these estimates are qualitatively consistent with 
the lattice results of Ref . , it should be added that they are dependent on the range of 
r-values used. 

This paper should be considered only as the first step in measuring densities. Many 
extensions and refinements are possible: 



In the Qg-system, the measurement of the Px/2, P3/2, D3/2, D$/2,--- densities 
corresponding to the energies extracted in Ref. 0. For a given orbital angular 



momentum, do these correlations show the degeneracy predicted in Ref. 

Off-axis radial correlations. These would check not only rotational invariance 
but also the radial correlation for r/a = 5, which could then be measured at, 
for example, x/a = 3, y/a = 4 etc. This point is of particular interest, since 
in Fig. ^| there is a hint that the correlation is lower than that obtained from 
a simple extrapolation. Such a lowering is expected, when the linear confining 
potential begins to play a role. 

The measurement of correlations in the baryonic and {Q 2 q 2 ) system. Are these 
similar to those in the (Qq) case - as is the case when comparing correlations in 
few-nucleon systems? 

The replacement of the quenched by unquenched gauge configurations. This is 
not expected to have a significant effect on the charge sum-rule and correlations. 
However, as discussed in Refs. [0], ||, for the matter probe disconnected contri- 
butions arise that are dependent on quenched versus unquenched. The difference 
between the two could highlight the effect of the quark condensate. 

The use of larger (3 values and lattices to check the continuum limit of the present 
results. 

The use of several light quark masses to enable an attempt at extrapolating to 
the realistic non-strange light quark masses. 

Other one body operators. In this work only the charge^) and matter (1) probes 
have been studied. However, others are of interest - see, for example, Ref. |T[ 
where the pseudo-vector operator (7^75) is needed for the B*Bn coupling. 
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APPENDIX A: THE CORRELATION FUNCTIONS C(2) AND C(3, R) 

This appendix gives details of the derivations of Eqs. |5] and |6] for the two- and three-point 
correlation functions C(2) and C(3). 

1. Two-point correlation functions C(2) 

Given the light-quark propagators in Eqs. and [5] , then the two-point correlation func- 
tions C(2) can be expressed in the four ways shown in Fig. [l]. These are the same upto 
statistical errors, but their combination improves the overall measurement. 

Case a: The Qq meson with the propagator in Eq. (Q). 

Case b: The Qq meson with the propagator in Eq. (|3|). 

Case c: The Qq meson with the propagator in Eq. (||). 

Case d: The Qq meson with the propagator in Eq. (||). 

In cases a and b, the heavy(static)-quark propagator from site (x, t) to site (x', t + T) 
is simply 

G (x, t ; x', t + T) = i(l + 7 4)t/ Q (x,t,T)5 x , x ,, (Al) 

where t/ Q (x, t, T) = Yii=o U±{x.,t + i) is the gauge link product in the time direction. On 
the other hand, for cases c and d, the heavy(static)-antiquark propagator from site (x, t) to 
site (x', t + T) is simply 

Gq(x, t ;x', t + T) = i(l- 74 )f/ Q t (X)t)T)5xx/ . (A2 ) 

The general form of a two-point correlation is constructed from a heavy-quark propagating 
from site (x 1; t) to site (x' 1; t + T) and a light-quark propagating from site (x 2 , t + T) to 
site (x' 2 , t) i.e. 

C(2,T) =Tr(r t G Q (x, £ ;x', t + T)rG,(x, t + T;x', t)), (A3) 

where T is the spin structure of the heavy-quark light-quark vertices at t and t + T. In this 
case, since the B- meson is a pseudoscalar, we have T = 75. These definitions lead to the 
four two-point correlation functions 

C(2, T,a)= £ (C/«(x, t, T)^(x, * + T)0,(x, *)) (A4) 

d=3,4 

C(2,r,6) = £ (f/ Q (x,i,T)0*(x,t + T)^(x,i)) (A5) 

(2=1,2 

C7(2, T, c) = J2 ( U * Q (*> *> T )^(x, t)0d(x, t + T)) (A6) 

(2=1,2 

C(2, T, d) = £ ( U * Q (*i t, T)02(x, tty d (x, t + T)), (A7) 

(2=3,4 
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where d is the Dirac spin index. We see that C(2,T,c)* = C(2,T,b) and C(2,T,d)* = 
C(2, T, a), so that the sum leads to 



C(2,T) = 2(Re 



U Q [J2 V:(x,t + T)0 e (x,t)+ £ 0^(x,t + T)^(x,t)] 



e=3,4 



d=l,2 



(A8) 



2. Three-point correlation functions C(3) 

The construction of the three-point correlation functions C(3) follows that of C(2) and 
are depicted in Fig. ||. 

Consider the probe is at t = and that Q propagates from (x, — £ 2 ) to (x, ti). Analogous 
to Eq. (|A3|) the general form of C (3) - when involving a Q - is then 

C Q (3,-t 2 , r)=Tr(r t G Q (x, -t 2 ;x, t^r G,(x, *i;x+r, 0) 6(r) Gj(x + r, 0;x, -t 2 )>. 

(A9) 

This combination of the G g and G' defined in Eqs. (|2]||) ensures that only the local field 
occurs at the probe position r. When 0(r) = 74, this expression reduces to 



C Q (3,-t 2 , t u r) 



d=l,2 



^ e (x + r, 0)<^* +2 (x + r, O)-0 e+2 (x + r, 0)<£(x + r, 0) 



e=l,2 



(A10) 



Care must be exercised here, since e (x + r, 0) arises from G q , whereas </4*(x + r, 0) is from 
G' q and so the two cannot be combined. The last equation can now be written as 



G«(3, -h, t u r) = ( [X q R + [if + iY?] ), 
where - supressing color indices - 



(All) 



e=l,2 



Xj = £ [0e(/)0l +2 (i?) " <f> e (R)<f>'e+2(I) ± <Pe + 2(R)<Pe(I) T ^W^)], ( A12 ) 



e=l,2 



where the upper signs are for O(r) = 74 and the lower ones for O(r) = 1. Supressing color 
and spin indices - 

Y§ = ^(R)U(R)^(R) + ^'(I)U(R^(I) - ^'(I)U(I)i;(R) + ^\R)U{I)^{I) 



Yj Q = -iP'(R)U(R)iP(I) + i[>'(I)U(R)i[>(R) + i[>'(I)U(I)i[>(I) + iP\R)U(I)^(R). (A13) 
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When involving a Q, the corresponding three-point correlation function is 

C§{-t 2 , h, r) = Tr(rG (x, t x ;x, -t 2 pG q (x, -t 2 ; x + r, 0) 0(r) G' q (x + r, 0,x; i x )>, 

(A14) 



which can be written in the notation of Eq. ( All 



as 



(A15) 



C«(3, -t 2 , t h r) = ( [-X% + iXfi [Y® - iY?\ ). 
The overall three-point correlation function is then 

C(3, -t 2 , t u r) = C Q (3) =F C Q (3) = 2( - X'Y? ), (A16) 

where the relative sign enters for 0(r) = 74, since the q and g have opposite charges. 
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TABLES 



TABLE I. Values of the parameters am a and vf, where a ps 0.17fm is the lattice spacing. Cases 
1 - 4 fit the two point correlations [C(2)]. 



am a 


Case 1 


Case 2 


Case 3 


Case 4 


vf uf 


3 3 T^rain — 3 


2 2 Train — 3 


3 3 Train — 4 


2 2 Tmin — 4 


ami 


0.8849(10) 


0.9005(14) 


0.8721(19) 


0.8833(27) 


am,2 


1.2953(63) 


1.355(10) 


1.263(13) 


1.307(20) 




1.99(10) 


— 


1.94(30) 


— 


< 


0.5164(30) 


0.5574(41) 


0.4847(56) 


0.5149(82) 


V F1 


1.5892(48) 


1.6761(52) 


1.519(10) 


1.589(13) 


V F2 


0.8651(22) 




0.8402(38) 




vl 

o 

V F 1 


0.8123(61) 
0.435(22) 


-0.8226(74) 
-0.065(29) 


0.816(16) 
0.644(49) 


-0.834(19) 
-0.373(57) 


vh 


-0.393(18) 


— 


-0.251(33) 


— 


vl 


-0.258(63) 


— 


-0.28(22) 


— 




1.93(32) 


— 


2.2(1.4) 


— 


vh 


-1.40(24) 


— 


-1.13(81) 


— 


ndatail) 


60 


30 


54 


27 


Tlparam (2) 

ndof(2) 
X 2 /n dof (2) 


12 

48 


6 

24 


12 
42 


6 

21 


3.1 


7.5 


0.65 


1.15 


< 


0.0421 


-0.0480 


-0.0608 


-0.3294 


U F1 


0.3546 


0.6126 


0.3258 


0.7361 


U F2 


0.4793 




0.6364 




i 


1.1135 


-1.2482 


1.1575 


-1.4024 


U F1 


-0.1432 


0.4151 


-0.1025 


0.4545 


U F2 


-0.4016 




-0.4825 
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TABLE II. Values of the x al3 for Cases 3 and 3'. In Case 3 the 12 parameters describing (7(2) 
are fixed and only the 6 x a/3 are varied to fit C(3). In Case 3' only the 3 x aa are varied. The 
superfix * implies that the number fixed at this value 







Case 3 






Case 3' 




X \2)/n dof (2) 




0.65* 






0.65* 






4 


6 


8 


4 


6 


8 


X 2 (3)/n dof (3) 


1.20 


0.72 


0.26 


4.2 


2.4 


0.68 


X 11 


1.26(2) 


1.33(2) 


1.41(5) 


1.12(1) 


1.19(1) 


1.29(3) 


X 12 


-0.32(3) 


-0.44(7) 


-0.5(2) 


0* 


0* 


0* 


x 13 


0.03(9) 


0.0(3) 


0(1) 


0* 


0* 


0* 


x 22 


0.65(7) 


0.9(2) 


0.9(9) 


0.00(4) 


-0.3(1) 


-1.2(4) 


x 23 


0.2(2) 


0(1) 


0(8) 


0* 


0* 


0* 


x 33 


0.2(9) 


4(8) 




-0.3(8) 


1(7) 








Case 4 






Case 4' 




X 2 (2)/n dof (2) 




1.15* 






1.15* 




Tmin(3) 


4 


6 


8 


4 


6 


8 


X 2 (3)/n do/ (3) 


1.23 


0.78 


0.21 


3.4 


2.0 


0.58 


x 11 


1.24(2) 


1.29(3) 


1.36(7) 


1.10(1) 


1.17(2) 


1.26(3) 


X 12 


-0.36(5) 


-0.45(11) 


-0.5(3) 


0* 


0* 


0* 


x 22 


0.65(13) 


0.8(4) 


1(2) 


-0.19(5) 


-0.6(2) 


-1.6(6) 
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TABLE III. The ground and first excited state charge densities x , x 22 for different values of 
T m j n (3). In Case 3 only the six x a ^(r) are varied - the v- vectors being those determined earlier 
from fitting C(2) and shown in Table |. In Case 3' only the three x aa (r) are varied. The numbers 
in [...] are the x 2 (3)/^do/(3). The entries denoted by '-' have unreasonably large x 2 (3)/^do/(3). 
The interquark distance r is in lattice units of a ~ 0.1 7fm. 



x^^ ,Case 


3 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


^min(3) 


6 


0.0627(2) 


0.0214(1) 


0.00906(8) 


0.00328(6) 


0.00112(6) 


0.00031(6) 






[0.37] 


[2.54] 


[2.37] 


[0.20] 


[0.40] 


[0.82] 


^mi?i(3) — 


7 


0.0629(3) 


0.0209(2) 


0.00882(10) 


0.00326(9) 


0.00120(9) 


0.00021(8) 






[0.27] 


[1.67] 


[1.49] 


[0.25] 


[0.34] 


[0.59] 


^min(3) = 


8 


0.0633(5) 


0.0204(2) 


0.0084(2) 


0.0033(1) 


0.0012(1) 


0.00017(11) 






[0.16] 


[1.38] 


[0.36] 


[0.33] 


[0.49] 


[0.68] 


^min(3) 


9 


0.0631(7) 


0.0200(4) 


0.0081(2) 


0.0033(2) 


0.0010(2) 


-0.0002(2) 






[0.24] 


[1.48] 


[0.05] 


[0.48] 


[0.53] 


[0.39] 


1 1 „ 

x , Case 


o/ 

3 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


-^min(3) — 


7 






0.00856(5) 


0.00287(4) 


0.00100(5) 


0.00018(5) 










[2.09] 

L J 


[1.87] 

L J 


[0.94] 

L J 


[0.58] 


-^min(3) — 


8 






0.0084(1) 


0.0029(1) 


0.0010(1) 


0.00012(7) 










[0.38] 


[1.08] 


[0.67] 


[0.81] 


Tmin(3) — 


9 






0.0082(1) 


0.0029(1) 


0.0009(1) 


0.0000(1) 










[0.07] 


[0.80] 


[0.37] 


[0.33] 


x 22 , Case 


3 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


T m i n (3>) — 


6 


0.192(3) 


0.011(1) 


0.0019(7) 


0.0009(5) 


0.0004(6) 


0.0003(6) 


^min(3) — 


8 


0.19(1) 


0.010(4) 


0.001(3) 


0.001(2) 


0.000(3) 


0.001(2) 


x 22 , Case 


3' 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


2~min(3) — 


7 






-0.001(1) 


-0.0033(4) 


-0.0020(6) 


0.000(1) 
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TABLE IV. The parametrization of the wavefunction as w c (r) = A c exp[— (r/r c ) p ] for various 
values of T m j n (3) and ranges of r. Also given is I = irA 2 r 3 . Here only fits with x 2 /dof ~ 1. The 
cases with p = 1* have p fixed. Both A c and r c are expressed in terms of lattice units a ~ 0.17fm. 



Case 3 


r— range 


A c 


r c 


P 


I 


^rmn(3) — 7 


2-4 


0.254(7) 


2.01(5) 


1* 


1.52(14) 




1-4 


0.199(5) 


2.87(14) 


1.23(5) 


1.52 


T m i n (3) = 8 


1-4 


0.225(3) 


2.20(4) 


1* 


1.71(10) 




2-4 


0.238(10) 


2.10(8) 


1* 


1.65(23) 




1-4 


0.208(11) 


2.59(25) 


1.13(9) 


1.56 


T m i n {2>) — 9 


1-4 


0.224(4) 


2.18(6) 


2 + 


1.64(15) 




1-4 


0.202(15) 


2.66(37) 


1.16(14) 


1.47 




1-5 


0.196(12) 


2.81(31) 


1.22(12) 


1.41 


Case 3' 


-^min(3) — 7 


2-4 


0.275(4) 


1.84(2) 




1.47(6) 




1-4 


0.249(5) 


2.12(8) 


1.08(3) 


1.43 


^min(3) = 8 


1-4 


0.257(2) 


1.94(2) 


2 + 


1.50(5) 




2-4 


0.261(6) 


1.91(4) 


1* 


1.49(12) 




1-4 


0.250(9) 


2.04(13) 


1.03(4) 


1.47 


^min(3) — 9 


1-4 


0.251(2) 


1.96(3) 


1* 


1.49(7) 




2-4 


0.260(10) 


1.90(7) 


1* 


1.46(20) 




1-4 


0.237(12) 


2.19(21) 


1.08(7) 


1.41 



TABLE V. Values of the x a @ in Cases 3 and 3' for the matter distribution. Notation as in 



Table || 







Case 3 






Case 3' 




X 2 (2)/n dof {2) 




0.65* 






0.65* 




-^min(3) 


4 


6 


8 


4 


6 


8 


X 2 (3)/n do/ (3) 


0.40 


0.45 


0.35 


0.72 


0.43 


0.29 


X 11 


0.56(3) 


0.49(6) 


0.38(15) 


0.48(2) 


0.45(3) 


0.35(7) 


X 12 


-0.22(7) 


-0.1(2) 


-0.1(5) 


0* 


0* 


0* 


X 13 


-0.2(2) 


0.0(8) 


-1(5) 


0* 


0* 


0* 


X 22 


0.4(2) 


0.4(5) 


1(2) 


-0.11(9) 


0.0(3) 


0(1) 


x 23 


0.3(4) 


0(3) 




0* 


0* 


0* 


X 33 


2(2) 






2(2) 
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TABLE VI. The Matter densities with notation as in Table 



III 



x 11 ,Case 3 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


^~min(3) = 6 


0.0758(3) 


0.0199(1) 


0.0062(1) 


0.00145(8) 


0.00040(8) 


0.00004(7) 




[0.20] 


[1.07] 


[0.83 ] 


[0.83] 


[0.23] 


[0.60] 


T m i n {3) = 8 


0.0758(7) 


0.0192(3) 


0.0054(3) 


0.0009(3) 


0.0004(2) 


0.0005(2) 




[0.25] 


[0.04] 


[0.19] 


[0.94] 


[0.19] 


[0.06] 


x ,Case 3 


r = 


r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


^min(3) = 6 








0.00113(4) 


0.00034 (4) 


0.00003(5) 










[1.82] 


[0.24] 


[0.57] 


T m i n {2>) = 8 




0.0196(2) 


0.0046(2) 


0.0010(1) 


0.00043(9) 


0.0003(1) 






[1.24] 


[0.80] 


[0.77] 


[0.16] 


[0.16] 


TABLE VII. The parametrization of the matter distribution wavefunction as 




w m (r) = A n 


exp[-r/r m ] - 


notation as in Table IV 








Case 3 


r 


range 


A 




' m 




T m i n {3) = 6 




2-4 


0.33(2) 




1.40(5) 


0.92(14) 


T m i n {2>) = 8 




0-4 


0.275(1) 




1.46(2) 


u. / o\o) 






1-4 


0.271(8) 




1.49(5) 


0.77(9) 






2-4 


0.37(8) 




1.24(16) 


0.81(46) 


Case 3' 


r 


range 


A 












1-4 


0.293(2) 




1.39(1) 


0.72(2) 


Tmin (3) = 8 




1-4 


0.287(5) 




1.39(3) 


0.69(5) 
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FIGURES 



FIG. 1. The four contributions to the two-point correlation function C(2) 




FIG. 3. a) The basic LL, F1F1 and F2F2 elements of C(2): b) The values of E(T) for LL, 
F1F1 and F2F2 separately: c) The combinations of L, Fl, F2 for cases 3(solid) and 4(dashed) to 
give E[Caa]. Also shown is the ratio R = {C X2 {2,T)/C l i{2,T) ) defined in Eq. || 




FIG. 4. The 74 sum rule: a) Contributions LL, F1F1 and F1F1 separately: b) The combina- 
tions of L, Fl, F2 for cases 3(solid) and 4(dashed) to give F sum [C aa (^4, T)] defined in Eq. 15 and 
R sum in Eq ija 
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FIG. 5. The ratio (C(3, r)/C(2)) for Case 3. a) r = 0,1,2 and b) r = 4,5,6 in lattice units 
a ~ 0.17fm. 
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FIG. 6. The radial distribution of the ground state charge (C) and matter(M) densities x n (r). 
a) Case 3(solid) and 3'(dashed) with T min (3) = 8: b) Case 3 for T m j n (3)=8(solid) and 6(dashed). 
The interquark distance (r) is in lattice units a « 0.17fm. 
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FIG. 7. Matter Sum rule with Cases 3(solid) and 4(dashed) L,F1,F2 combinations 
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